Self-similarities in the frequency-amplitude space of a loss- modulated CO2 laser 
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We show the standard two-level continuous-time model of loss-modulated CO2 lasers to display 
the same regular network of self-similar stability islands known so far to be typically present only 
in discrete-time models based on mappings. For class B laser models our results suggest that, more 
than just convenient surrogates, discrete mappings in fact could be isomorphic to continuous flows. 
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Lasers with modulated parameters are arguably among 
the simplest and most accessible laser systems of inter- 
est for applications in science and engineering and for 
theoretical investigations. The intrinsic interest in prac- 
tical applications and in the nonlinear dynamics of mod- 
ulated lasers has spurred a wide range of studies after the 
remarkably influential work of Arecchi et al. report- 
ing the first measurement of subharmonic bifurcations, 
multistability, and chaotic behavior in a Q-switched CO2 
laser. Since then, CO2 lasers have been fruitfully ex- 
ploited in many situations. Recent applications include 
studies of stochastic bifurcations in modulated CO2 laser 
0, multistability induced by periodic modulations 
Rich nonlinear response of CO2 lasers with current mod- 
ulation and cavity detuning and self- focusing effects 
in nematic liquid crystals 

In the last 20 years the CO2 laser was extensively 
studied theoretically, numerically and experimentally 

0,11, 01 but focusing mainly on the characterization of 
dynamical behaviors in phase-space for specific parame- 
ters. While a detailed description of phase-space dynam- 
ics is already available in the literature H, 0, 0, 0| , 
no equivalent description exists for the parameter space, 
except for works by Goswami [l^ who investigated ana- 
lytically the first few period-doubling bifurcations for the 
Toda model of the CO2 laser [T3 |. 

The present Letter reports an investigation of the pa- 
rameter space of a paradigmatic model of class B lasers, 
the CO2 laser. More specifically, we study a popular 
two-level model of a CO2 laser with modulated losses, 
focusing on the global stability of the laser with respect 
to the modulation, not the intensity. The remarkable 
discovery reported here is that stability islands of the 
continuous-tmic laser model emerge organized in a very 
regular network of self-similar structures called shrimps 
|15j . illustrated in Figs. ^ and [SJ and previously known 
to exist only in the par ameter space of discrete-tiuie dy- 
namical systems [la [l^ 113, llal • Thus far, all attempts to 
uncover shrimps in flows, i.e. in continuous-time dynam- 
ical systems modeled with sets of differential equations, 
have failed to produce them [T^ . 

The single-mode dynamics of the loss-modulated CO2 
laser involves two coupled degrees of freedom and a time- 



dependent parameter which we write, as usual 
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Here, u is proportional to the radiation density, z and 
zo are the gain and unsaturated gain in the medium, re- 
spectively, T denotes the transit time of the light in the 
laser cavity, 7 is the gain decay rate, and k = k{t) rep- 
resents the total cavity losses. The losses are modulated 
periodically as follows. 



k{t) = ko{l + acos27r/i). 



(2) 



where fcg is the constant part of the losses and a and 
/, the amplitude and frequency of the modulation, are 
the main bifurcation parameters. The remaining param- 
eters are fixed at r = 3.5 x 10"^ s, 7 1.978 x lO'^ 
s~^, zq = 0.175 and kg = 0.1731. These are realistic val- 
ues, used in recent theoretical and experimental investi- 
gations Q- Integrations were done using the standard 
fourth-order Runge-Kutta scheme with fixed time-step, 
equal to /i = 4 X 10^^. Phase diagrams in a x / space are 
obtained by computing Lyapunov exponents for a mesh 
of 600 x 600 equally spaced parameters. Starting from 
an arbitrary initial condition, we "followed the attractor" 
that is, after increasing parameters we initiated iterations 
using the last obtained values as the new initial condi- 
tions. The largest exponents were codified into a bitmap 
with a continuous color scale ranging from the maximum 
positive (green) to maximum negative (blue) exponents. 
Zero exponents were codified in black. One of the expo- 
nents is always zero since it is simply related to the time 
evolution. Three illustrative bitmaps for the laser model 
are shown in Fig. ^ 

Figure^ displays a global view of the parameter space. 
The most prominent features, the broad curved struc- 
tures in Fig. show that the parameter space of the 
laser model above, Eqs. l^-b), agrees qualitatively quite 
well with the description of Goswami [l3| for the Toda 
model of the CO2 laser. For the parameters chosen, the 
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FIG. 1: (Color online) Structure of the frequency-amplitude phase-diagram of the laser showing a regular network of stability 
islands (a) Global view, (b) Zoom of the box in (a). Numbers indicate the main period of each stability island; (c) Magnification 
of the period-11 stability island (indicated in (b) by the encircled 11), displaying the generic shape of all stability islands [l^ . 
Color intensities are proportional to Lyapunov exponents: blue for negative exponents (periodic oscillations), black for zero 
and green for positive exponents (chaotic oscillations). Frequencies are in kHz. 




FIG. 2: (Color online) Structure of the parameter space of the Henon map. (a) The organization of shrimps here coincides 
with that of the laser (compare with Fig. 03). The fine structure observed around period-8 here reproduces the laser period-8 
sequence along the curve passing by the encircled numbers in Fig. 0d; (b) Magnification of the box in (a); (c) Magnification 
of the box in (b). Numbers indicate the main period of each stability island. Points mark windows, not doubly superstable 
crossings. Download ALL figures high-res PDF: |http:/7 www . if . ufrgs . br / ~j gallas/ j g_papers . ht ml j 



relaxation frequency of our laser model is 50 kHz. From 
Fig. Hi it is possible to see that there is a minimum ampli- 
tude threshold a beyond which subharmonic bifurcations 
start to occur, corresponding to about 100 kHz, the har- 
monic of the relaxation frequency. In addition, for certain 
parameter values new stability domains are created by 
saddle-node bifurcations, each of them undergoing then 
its own cascade of period doublings. So, in certain pa- 
rameter ranges more than one stable mode coexist, giv- 
ing rise to generalized multistability. This feature may 
be recognized in Fig. ^ from the apparent sudden discon- 
tinuities in the coloring, due only to the impossibility of 
plotting two distinct colors in the same place. 

The most interesting feature in Fig. ^ is the remark- 
ably regular structuring which appears in the region con- 
taining the box, shown magnified in Fig. ^p. This fig- 
ure shows that embedded in the wide domain of parame- 



ters leading to chaotic laser oscillations there is a regular 
structuring of self-similar parameter windows, shrimps, 
containing cascades of stable periodic oscillations, the 
main period of a few of the larger shrimps indicated by 
the number near to them. The period-11 shrimp seen 
in Fig. is shown magnified in Fig. Starting from 
the main period-11 body, it displays two distinct dou- 
bling cascades as well as an infinite number of addi- 
tional period-doubling cascades, as thoroughly described 
for discrete-time systems in Refs. [TslIT^. 

The computation of bitmaps for the laser model is very 
computer demanding. To alleviate this problem and to 
manifest the isomorphism between fiows and maps, we 
display the generic fine and hyperfine structure of stabil- 
ity islands typically present in multidimensional systems 
using the two-parameter Henon map as a paradigm: 

xt+i = a- Xf +byt, yt+i = Xt- (3) 
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FIG. 3: Comparison of return maps. Left column: Laser 
return maps for period-8 and its pair of doublings. Parame- 
ters are: (a) (a,/) = (0.06984, 89.8), period 8, (b) (0.07138, 
90.47), period 16, (c) (0.06902, 87.43), period 16. Frequencies 
are in kHz. Right column: Henon return maps for period-8 
and its doublings seen at the center of Fig. Parameters 
are: (d) {a,b) = (1.8087, -0.02514), period 8, (e) (1.80395, 
-0.0257), period 16, (f) (1.80642, -0.02356), period 16. 
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FIG. 4: Time evolution of the laser intensity u{t) for the large 
period-8 structure in Fig.^, plotted in (a) logarithmic scale, 
(b) linear scale. Here T = 1/(89.8 kHz) is the period of the 
modulation. 



The nonlinearity parameter a (forcing) represents the bi- 
furcation or control parameter. The damping parameter 
b varies between — 1 < & < 1 , with 6=1 representing the 
conservative limit and 6 = the limit of strong damping. 
While for 6 = there exists just a single chaotic attrac- 
tor over a wide range of the parameter a, several periodic 
and chaotic attractors coexist when 6 7^ 0. 
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FIG. 5: Predicted signal intensity for the six period-16 sta- 
bility islands labeled A, B, . . .F, in Fig. Signals (a) and 

(b) , at points A and B, are period-8 doublings. All other sig- 
nals are from islands which begin with period-16. Note that 
signals look very similar, despite the fact that they originate 
from very distinct regions of the parameter space. Parameters 
(a,/) are: (a) A = (0.06902, 87.43), (b) B = (0.07138, 90.47), 

(c) C = (0.063715, 92.15), (d) D = (0.062255, 75.735), (e) 
E = (0.0749617, 67.3281), (f) F = (0.073666, 83.359). Note 
that T = 1// is slightly different for each signal. 



Since all these attractors evolve in the same generic 
manner, we consider here the strongly dissipative limit, 
focusing on slightly negative values of 6. In this domain, 
Pando et al. |2(J | found that a much more sophisticated 
four-level model of the CO2 laser with modulated losses 
behaves qualitatively similar to the Henon map. The 
CO2 laser dynamics, as of any class B laser, is charac- 
terized by a time-delay between the intensity and the 
population inversion, a fact that nicely matches the de- 
layed character of the Henon map when written as a one- 
dimensional recurrence relation. 

Figure 121 shows how regularly shrimps organize them- 
selves along very specific directions in parameter space. 
The ordering along the main diagonal of Fig. is the 
same found for the laser, in Fig. along the direction 
containing the encircled periods. Analogously, the sec- 
ondary diagonal in Fig. |2K displays the same ordering 
that the parabolic arc in the middle of Fig. 

The laser-Hcnon agreement in parameter space perme- 
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ates also to the phase space as corroborated by Fig. |31 
comparing return maps for the laser (on the left column) 
with those for the Hcnon map (right column). The laser 
return maps were constructed using the sequence ug{t), 
£ = 1,2,3, .. . of relative maxima of u{t). As it is easy to 
see, both sets of return maps agree remarkably well |2ll |. 

How easy is to detect experimentally the regular struc- 
turing reported here? Figure ^ illustrates a represen- 
tative laser signal in two scales. Although waveforms 
and underlying periodicities are easy to recognize in log- 
arithmic scale, their experimental detection may become 
strenuous, particularly as the period increases. For in- 
stance, contemplating the six period-16 stability islands 
in Fig. njj, two of them arising from period-8 via saddle- 
node bifurcations, one may ask what sort of differences 
should be expected in their measurement. The answer 
is depicted in Fig. |S| In a real experimental setup, the 
difficulties to surmount are mainly to access narrow high- 
period windows, and to have a wide enough detection 
range. Modulated losses are usually obtained with an 
intracavity polarizer and an electro-optical modulator. 
Recent progress in low- voltage electro-optical modulators 
have considerably improved their performances p3 |. To 
detect large and small peaks simultaneously one can use 
a logarithmic preamplifier Thus, detection and dis- 
crimination of the laser signals in Fig.jSjis experimentally 
feasible with existing technology. 



To uncover isomorphisms between continuous-time 
(flows) and discrete-time (maps) in dynamical systems is 
an important event from a physical and dynamical point 
of view and immediately raises interesting questions. For 
instance, is an isomorphism to be expected also for more 
refined laser models such that, e.g., of Ciofini et al. [23| . 
involving two rate-equations derived for a single-mode 
CO2 laser using center manifold theory for a four-levels 
model (2^? Their model is interesting because, as they 
say, it agrees well with experiments. Do networks of 
shrimps exist for the laser model proposed very recently 
by Meucci et al. . Do they exist for other laser models 
too? Is it possible to find them in the parameter space 
of the laser parameters, not modulation? In autonomous 
systems? Which laser system presents more restricted 
amplitude variations, making life easier for experimental- 
ists? What sort of mechanism generates shrimps in differ- 
ential equations? How are basins of attraction entangled 
in multistable domains where self-similar structures are 
abundant |23|? We will report on this elsewhere. 
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